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Abstract. 

Using the IRAM interferometer we have mapped at high 
resolution (2'.' 2 x 1'.'2) the CO (1-0) emission in the nucleus 
of the doubled barred SABbc spiral M 100. Our synthesized 
map includes the zero spacing flux of the single-dish 30m map 
(Sempere & Garcfa-Burillo, 1997, paper I). Molecular gas is 
distributed in a two spiral arm structure starting from the end 
points of the nuclear bar (r = 600 pc) up to r = 1.2 kpc, and a 
central source (r ~100 pc). The kinematics of the gas indicates 
the existence of a steep rotation curve (vrot= 180 km at 
r ^ 100 pc) and strong streaming motions characteristic of a 
trailing spiral wave inside corotation. 

Interpretation of the CO observations and their relation 
with stellar and gaseous tracers (K, optical, Ha, HI and radio- 
continuum maps) are made in the light of a numerical model 
of the clouds hydrodynamics. Gas flow simulations analyse the 
gas response to a gravitational potential derived from the K- 
band plate, including the two nested bars. We develop two 
families of models: first, a single pattern speed solution shared 
by the outer bar+spiral and by the nuclear bar, and secondly, 
a two independent bars solution, where the nuclear bar is dy- 
namically decoupled and rotates faster than the primary bar. 

We found the best fit solution consisting of a fast pat- 
tern {^If — 160 kms^^kpc"^) for the nuclear bar (with coro- 
tation at Rqqji=1.2 kpc) decoupled from the slow pattern of 
the outer bar-l-spiral (f2/=23 kms~^kpc~^) (with corotation at 
R^Q^=8-9 kpc). As required by non-linear coupling of spirals 
(Tagger et al 1987), the corotation of the fast pattern falls in 
the ILR region of the slow pattern, allowing an efficient transfer 
of molecular gas towards the nuclear region. Solutions based 
on a single pattern hypothesis for the whole disk cannot fit the 
observed molecular gas response and fail to account for the 
relation between other stellar and gaseous tracers. In the two- 
bar solution, the gas morphology and kinematics are strongly 
varying in the rotating frame of the slow large-scale bar, and 
fit the data periodically during a short fraction (about 20%) 
of the relative nuclear bar period of 46 Myr. 
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1. Introduction 

The advent of high-sensitivity near-infrared imaging of galax- 
ies has shown that a significant percentage of barred spirals 
host secondary bars in their nuclei. There could be two inter- 
pretations of the bars within bars phenomenon, according to 
the relative pattern speeds of the two bars (Friedli and Mar- 
tinet, 1993; Friedfi and Benz, 1993 and 1995; Combes, 1994). 
The two patterns could be corotating if they are about parallel 
or perpendicular to each other. If the secondary inner bar is 
strongly misaligned with the primary outer bar, the two bars 
are likely to have distinct wave pattern speeds, as shown by 
numerical simulations. The decoupling of an inner faster pat- 
tern appears in self-consistent simulations with gas and stars 
thanks to the role of the dissipative component: as a result of 
gas inflow, under the action of the bar gravitational torques, 
mass accumulates onto the X2 types of orbits which weakens 
the primary bar. The rotation period becomes much shorter in 
the nuclear regions due to mass concentration, which leads to 
the decoupling of a fast-rotating bar. Eventually, the nuclear 
bar destroys itself or it destroys the primary bar, modifying 
the overall disk potential. Evolution can then occur in much 
less than a Hubble time, and galaxies change their morpholog- 
ical type along the Hubble sequence. They change from barred 
to un-barred, and also they concentrate mass in the process, 
evolving slowly from late-types to early types. 

The observation and modelling of real barred galaxies offers 
the opportunity to test theory predictions on galaxy evolution 
and it seems a necessary complement to numerical simulations 
of model galaxies. The present work is intended to bring a com- 
bined observational and modelling effort in the study of the 
nearby barred spiral MlOO (NGC4321). In this galaxy, classi- 
fied as SABbc by de Vaucouleurs et al (1991), the hypothesis 
of a single mode common to the whole disk is dubious, both 
from observational and theoretical evidences. 

On the observational side, the nuclear region of MlOO (up 
to r=3kpc) has been so far the subject of numerous stud- 
ies. The pioneering work of Arsenault and collaborators (1988, 
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1989, 1990) established a connection between the ring-like Ha 
morphology of the nucleus and the existence of Inner Linblad 
Resonances. Further steps in sensitivity made appear, first, a 
four-armed structure (Cepa et al, 1990) and recently a frag- 
mented two spiral arm structure in Ha (Knapen et al, 1996). 
The 6cm radio-continuum VLA map of Weiler et al (1981) 
shows also a two arm spiral pattern. 

Near infrared images of the nucleus (Pierce 1986, Shaw et al 
1995, Knapen et al 1995 (hereafter K95), Rauscher 1995) show 
the existence of either a secondary nuclear bar in K (nearly par- 
allel to the main bar) together with a leading spiral structure, 
or an inner oval in the I band (with principal axes misaligned 
with respect to the outer bar). The synthesis aperture ^^CO(l- 
0) maps (Canzian, 1992; Rand, 1995; Sakamoto et al 1995) 
indicate the existence of a two-arm molecular spiral structure 
connected to the K nuclear bar end points. The IRAM 30m 
map of paper I shows a strong concentration of CO emission 
towards the nuclear disk ND, a component clearly distinguish- 
able from the main bar. A steep rotation curve gradient, unre- 
solved by the 30m beam (12" in the 2-1 line), indicates a high 
mass concentration in the ND. 

On the modelling side, Garcia-Burillo et al (1994) (here- 
after called GB94) and Sempere et al (1995) (hereafter S95) 
made numerical simulations of the cloud hydrodynamics to 
study the evolution of the molecular gas disk under the ac- 
tion of a realistic spiral+barred potential derived from a red 
band plate. The authors assume the whole disk to be fitted by 
a single well defined wave pattern characterized by Qp, shared 
by the primary bar and the spiral arms. However they lacked 
first, of the necessary spatial resolution and secondly, of a fair 
potential tracer to analyse the gas response in the inner 500 
pc. 

K95 have made numerical simulations of the stellar and 
gas dynamics in MlOO, using a model potential which de- 
parts markedly from the real mass distribution. Although they 
favour a one bar mode scenario their model fails to reproduce 
the molecular gas distribution observed by the interferometer. 

We present here a combined single-dish and interferometer 
data set fulfilling both high resolution (2.2" x 1.2") and sensi- 
tivity requirements. Contrary to the synthesis aperture maps 
so far published, we recover entirely the zero-spacing fiux of the 
ND. The comparison between the different gaseous and stellar 
tracers of the ND is reexamined in this work. Particular at- 
tention is paid to the bias introduced by extinction in optical 
and even near-infrared images of the nucleus, and what might 
be the implications on the interpretation of the data. Observa- 
tions are confronted to the result of new numerical simulations 
of the clouds hydrodynamics, based on a mass distribution di- 
rectly derived from the infrared luminosity image of MlOO. We 
will focus on the feasibility of two independent patterns in the 
disk and how this scenario accounts better for the observations. 

2. Observations 

^^CO(l-O) visibilities were obtained between February 1995 
and April 1996 by the IRAM Plateau de Bure interferometer. 
We used the 4-antennas set of extended configurations (BC) 
supplemented by 3-antennas observations compact set (D). All 
antennas were equipped with SIS receivers operating in the 
single-side-band mode and yielding receiver temperatures close 



to 40 K. The system temperatures range typically from 250 K 
to 400 K. The cross-correlator was adjusted within a 540 MHz 
passband with a channel width of 2.5 MHz, corresponding to a 
velocity resolution of 6.5kms~^ in the ^^CO(l-O) line. 

Instrumental calibration was performed using the quasars 
3C273 and 1219-1-285. The atmospheric rms phase fluctuations 
on the longest baselines, integrated over 4 minutes, were al- 
ways found to be < 50 degrees. The receiver passband was 
verified on 3C273. The absolute flux density scale was based 
on measurements of 3C273 and 1219-1-285. Information on 
their fluxes were bootstrapped from regular IRAM monitor- 
ing observations (Dutrey & Ungerechts, 1995, 1996). The data 
were antennar-based calibrated using the CLIC package (Lucas, 
1992). 

Short spacing visibility data was provided by the 30m 
radiotelescope. Single-dish visibilities were derived from the 
brightness temperature distribution following a reduction 
scheme described in the software package GILD AS (Guilloteau 
& Forveille, 1989) (see also the comprehensive discussions in 
Vogel et al. 1984, Neri et al. 1997). Firstly, the single-dish data 
were rcsampled to the 2.5 MHz frequency resolution of the in- 
terferometric data, and the main beam temperatures were con- 
verted to Jy/beam by applying the standard conversion factor 
S/TX = 5.5 Jy/K adopted for the 30m telescope (Wild 1995). 
Secondly, the single-dish maps were deconvolved by dividing 
their Fourier transforms by the Fourier transform of the 21" 
beam of the telescope, which was assumed to be circular and 
Gaussian, and the results multiplied by the 42" primary beam 
of the interferometer. Finally, the single-dish map was cor- 
rected from a 2" absolute residual pointing error. The derived 
MlOO dynamical center (sec below) coincides with the emis- 
sion peak in the single-dish and the interferometer synthesized 
field. 

Single-dish and interferometric visibilities were then com- 
bined into a common visibility table by scaling the single-dish 
weights to the mean weight of the short spacing interferometer 

visibilities in the uv range 15 — 25m. Henceforth we implicitly 
correct for eventual calibration errors in the 30m data. 

CLEANed maps were then obtained from the visibilities by 
standard deconvolution procedures (no tapering and uniform 
weighting). The map field is 256 x 256 pixel^ in size with a 
linear scale of O'.'25/pbcel. The synthesized beam was 2" 2 x 
l"2 and oriented along PA=30 degrees. The rms noise level 
in the 2.5 MHzchannel maps is ^ 4mJy/beam or oquivalently 
140mK within the synthesized beam. Point source sensitivity, 
determined for the line-free 400 MHzcontinuum passband is ~ 
0.65 mjy/beam. 

We adopted a(1950)= 12''20'"23.2% <5(1950)= 16°06'00"as 

the phase tracking center. Correlator was centered at 
vls-R = 1562 kms^"'^. However, we re-dctcrniincd the dynamical 
center and the systemic velocity of NGC4321, following a stan- 
dard least squares procedure applied to the interferometer data 
cube. Assuming a priori a position angle of PA=153° (paper I), 
the dynamical center (the (0,0) offset as referred in this work) 
is found at a(1950)=12''20"'22.96^ 5(1950)=16°05'57.5"; the 
systemic velocity is fitted to Vsys{LSR) — 1573iiz3kins^^. The 
velocity and spatial scales used in this paper will be relative to 
the derived VsysiLSR) and (0,0) central offset. At the assumed 
distance of MlOO {D ~ 17.1 Mpc), 1" ~ 82 pc. 
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3. Molecular gas distribution 

3.1. The CO maps 

We present in Fig.l the ^^CO(l-O) velocity-channel maps dis- 
playing line emission from v=-91knis^^ to v=91kms^^. No 
significant emission is found outside this velocity range. The 
rms noise in the channel maps (6.5kms~^ wide) is 8.5mK in 
brightness temperature. The emitting gas is confined within 
an elongated region extending over Ax x Ay =20" x 30". As ex- 
pected for a moderately inclined (i=32°) rotating disk whose 
major axis is close to North (PA=153°) and assuming that 
MlOO's spiral structure is trailing, gas emission is bound to 
appear mostly at positive velocities on the southern side, con- 
trary to blueshifted gas emission showing up on the northern 
side. Notwithstanding this general trend, the morphology of 
the channel maps indicate strong deviations from axisymmetry 
both in the velocity field and in the gas distribution. Channel 
maps with velocities close to Ysya=Q, show a distorted S-like 
pattern reminding of spiral-like streaming. Emission from a 
strong CO source is detected within a large velocity range 
v=(-70 up to 70kms^^) at the center of the galaxy. There 
arc also evidences of gas emission at highly forbidden veloc- 
ities: CO emission is detected on the northern part, towards 
(Aa, A5)= (-|-7",-|-10"), at positive velocities (v=60 up to 
SOkms"'^). There is a southern counterpart of this component 
towards (Aa, A(5)= (-7",-10"), at negative velocities (v=-60 
up to -80kms~^). 



Fig. 1. i^GQd 
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de Bure array with a resolution (HPBW) of 2.2"x 1.3". The dy- 
namical center in (1950) equatorial coordinates is indicated by a 
cross at a= 12^ 20™ 22?96 , 5= 16° 05' 57.5". We display the emis- 
sion in 6.5kms~^-wide velocity channels from v=-91kms~^ up to 
v=91kms~^, symmetrically disposed with respect to the systemic 
velocity v=0 (vs!/s=1573kms~i(LSR)). Contour levels are -20, 20 
to 180 mjy/beam by steps of 20 mjy/beam. 



The CO velocity-integrated maps (not corrected for pri- 
mary beam attenuation) for the Plateau de Bure (Ico) and 
the combined Plateau de Bure-l-30m data {Vco^) are shown in 
Figs2a-b. They have been obtained by integrating the emission 
in channels with brightness temperature higher than 2a. The 
morphology of molecular gas distribution seen by the interfer- 
ometer is characterized by: 



■■„^^ig, 2. a (top) ^■^CO(l-O) velocity-integrated intensity con- 
tpuis observed within the central HPBW=43" interferometer field 
'oT I M IOO. Aa and A5 offsets (in arcsec) are relative to the 
dynamical center. Gray level contours correspond to 0.25, 0.5 
and IJykms^^bcam^^ to 10.5 Jykms~^beam~^ by steps of 
0.75 Jykms~^beam~^. Ellipses trace the loci of resonances in the 
disk (see text) b (bottom) Same as a, but for the combined Plateau 
de Bure-|-30m map. We indicate the positions of clumps I-IV and C 
(see text). 



Two spiral arms stretching out from radius r'^7"and po- 
sition angles PA=130°(on the eastern side, hereafter Arm 
E) and 310° (on the western side, hereafter Arm W). The 
arms spread over an azimuthal range of A$=120°and reach 
the edges of the map. Arm W has a sudden change in its 
pitch angle near PA=210°, appearing broken and discon- 
tinuous. CO emission along the spiral arms is uneven: it 
displays a clumpy hierarchy. Some of these clumps (denoted 
as I-IV) have been spatially resolved by the interferome- 
ter beam. Spiral arms are narrow: at places, they are not 
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resolved transversely (FWHM< 2.2"). The sign of the ob- 
served non-circular motions both in the intcrarm and in 
the arm regions correspond to the expected velocity pat- 
tern inside corotation (see discussion of S95 and section 
6). 

— There is a strong concentration of molecular gas at the 
center in the shape of an ellipsoidal source which is 
marginally resolved by the interferometer beam (major 
axes-(a,b)=(2.5", 2.1")). The central clump C is connected 
to the western spiral by a bridge of emission. 

— Away from the spiral-like and central sources, gas emission 
is hardly detected in the inner and outer interarm region. 

Fig 2b shows the Iffo'' map. Roughly the morphology of 
both majjs arc alike. The zero-spacing flux recovered (defined 
as (f=Ioo''-Ico)/Ico reaches on average ~50%. However, 
this percentage varies markedly in the map: it lowers up to 10- 
20% for Arm E, it rises up to 50% for Arm W, and it reaches 
100% for the interarm region. The flux recovered in the central 
source (f~30%) makes it appear rounder and spatially resolved. 
Spiral arms are also broader and there is smooth intcrarm emis- 
sion detected up to r=10"that the interferometer map filtered 
out completely. Still, the ND arm-interarm contrast is high 
(~4-5). 

3.2. Mass of molecular gas 

We derive masses and column densities of molecu- 
lar gas from Ico using a CO-to-H2 conversion factor 
X=N(H2)/Ico=2.3xl0^°cm-2K-ikm-is (Strong et al. 1988). 
Assuming the distance to be D=17.1Mpc, the total molecu- 
lar gas mass in the interferometer and combined integrated 
intensity maps are 1.4xlO^M0 and 2.7xlO®M0 respectively. 
The masses of clumps (I, II, III, IV and C: all of thcni of 
sizes spatially resolved by the interferometer beam) are esti- 
mated, first, from I^'q '^'i Ice?'' (CO based masses including 
the correction factor for helium, Mco) and, secondly, from the 
virial theorem (Mvir). In the latter case, masses come from 
M(Il2)=550dcri^, where d is the equivalent FWHM diameter, 
defined as d=(difi2)^''^ (di and d2 are the FWHM diameters 
along the principal axes in a 2D-gaussian model), and aid is 
the measured radial velocity dispersion seen in the CO line. 
Clumps I, II, III and IV seem gravitationally bound entities 
(Mco = l-7-2.7Mv/fl, according to mass estimates from Vao^), 
contrary to the central clump (C), as expected (for C the mea- 
sured aid is entirely due to the large velocity gradient of ro- 
tation curve). CO-based masses range from 0.6 to 2.4x10® M0 
typical of GMA-like associations. Note however that we have 
no accurate calibration of the conversion factor in Ml 00 and 
that the previous conclusions on gravitational boundness of 
GMAs lean on the assumed value of X. 

The I'cQ contours at half power along the spiral arms cor- 
respond on average to molecular gas column densities of 400- 
SOOAfo/pc^. N(H2+He) peaks at SOO-lOOOMo/pc^ along both 
spiral arms and it reaches 2OOOM0/pc^ towards C. The criti- 
cal column densities for gravitational instabilities (Nc; Toomre 
(1964), see also Kennicutt (1989) can be derived from (see pa- 
per I) 

Nc = aavK/G 

Nc range from 100 (at r~15") to SSOMq/pc^ (towards C), 
and henceforth the condition for the onset of gravitational in- 



stability (N(H2-I-He)/Nc >1) is largely fulfilled along the spiral 
arms and the rmclcar source (N(H2+He)/Nc ~2 6) 

We used the relation between Ii2CO(i-o) and At, obtained 
by Cernicharo and Guelin (1987) for the Taurus cloud complex 
to estimate the visual extinction towards the center of MlOO: 

I('^CO(1-0))=(5.0±0.5)(2A„-0.5±0.2) [1] 

where the factor of 2 corrects for the non-obscuring dust 
behind the HII regions, assumed to be located in the middle 
of the disk. A„ derived from [1] ranges from 2-2.5 (in the 
gaseous spiral arms) to 5 (towards C) in good agreement with 
the average value obtained from I("CO(1-0)) in paper I for 
the ND (A„ - 3.5). 

4. Comparison with other tracers 

We discuss in this section the global picture emerging from the 
comparison of the different gaseous and stellar tracers in the 
center of MlOO. 

The near infrared K pictures of the MlOO disk show the 
presence of an outer stellar bar with a 90-100"diameter. A 
gaseous molecular bar, aligned with the stellar bar, has been 
detected (paper I) and a flow of gas parallel to the bar is 
evident in the observed CO kinematics up to r~60": the CO 
p-v minor axis diagram shows the S-like pattern typical of a bar 
driven gas flow (GB94). The molecular gas bar has two ridges 
offset from the bar major axis, mimicking the offset leading 
dust-lanes often observed in barred galaxies. This indicates the 
presence of two ILRs, and the tendency of the gas flow to 
become perpendicular to the bar in between (Athanassoula 
1992). 

The previous infrared (I, J and R) images of the nuclear 
region showed the existence of a imclear oval or inner bar, 
but with a 30° leading angle respect to the outer bar. The 
discrepancy is explained by differential extinction, affecting 
less severely the K band image. The spiral-like distribution 
of molecular gas and henceforth the extinction maxima locus 
explain why Pierce (1986) see a misaligned nuclear lens. K95 
used the alignment of both K bars as an argument support- 
ing the existence of one mode for the disk, opposed to the two 
independent bars scenario. Still, this is just a statistical argu- 
ment: whatever relative orientation of two independent bars is 
equally probable, including the alignment case. On the other 
hand, although the K band picture of MlOO is less affected by 
extinction, the existence of substructures such as leading arms 
or relative maxima at the bar ends (used by Knapen and col- 
laborators as supporting evidences of the one mode case), are 
shown to be linked with extinction (see below). 

4.1. Bias of extinction 

Fig. 3a shows the overlay of Pqq contours with H„, in gray 
scale. Although, similarly to CO, Hq emission is confined in a 
ring-like structure and a central source, there are systematic 
shifts between the CO and ridges along the arms. Assum- 
ing the spiral to be trailing (see section 6), Ha appears up- 
stream the CO ridge along both arms for PA=0°-90° (Arm E) 
and PA=180°-270°(Arm W). On the contrary, the H^ ridge 
lies downstream the CO maxiirmm in the aziirmthal ranges 
PA=90°-180°(Arm E) and PA=270°-360° (Arm W). Towards 
C, He is relatively dimmed. Generally speaking, there is ap- 
parently a local anticorrelation between molecular gas column 
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densities and massive star formation in the ND: CO peaks 
suspiciously avoid He maxima. 

From the map we can estimate the opacity in the Ha 
band and correct the map of Fig.Sa by a factor e^^a . Although 
some of the relative shifts between CO and Ha clumps axe still 
present after correction, many of the reported displacements 
disappear or at least they arc considerably attenuated. The 
latter result indicates that a direct comparison between CO 
and extinction uncorrected H^ maps may be highly misleading, 
as expected, considering that A„ ranges from 2 to 5, as shown 
above. 

Additional insight into extinction effects is illustrated in 
Fig.Sb, showing the overlay between Y^q (gray scale) and ra- 
diocontinuum emission at 6cm (Seem) (contours). MlOO pos- 
sesses an extended source of radio emission associated with the 
ND. The brightest radio peak lies at a (1950)=12''20'"23.46^ 
(5(1950)= +16°05'56.5". A secondary radio maximum coincides 
with the optical nucleus. Seem contains both thermal free-free 
emission and non-thermal synchrotron emission. In order to 
estimate the contribution of each component we have consid- 
ered, first, the existent measurements at three wavelengths 
(1.2cm, 6cm and 20 cm) and, secondly, a simple hypothesis 
on the spectral indexes of thermal and non-thermal emission 
(at/i=-0.1 and Q„onth=-1-0, respectively, where a is defined as 
Su ~ f~"). The resulting best fit indicates that thermal emis- 
sion from HH regions dominates at 6 cm, as Sscm.tft, >60%. 
As shown in Fig.Sb, there is an excellent spatial correlation 
between If'o and Seem- Henceforth, no systematic ofi^sets exist 
betwoou massive star formation, traced by Seem, (not affected 
by extinction) and molecular gas column densities, traced by 

Tint 
'-CO- 

Fig. 4 shows the overlay of Vqq contours with the K-band 
imago (in gray scale) which represents the bulk of the old stel- 
lar component. The isophote deviations in the inner 6" were 
first interpreted by K95 as the signature of leading arms. A 
similar interpretation is invoked to account for the K peaks 
visible at the bar end-points (denoted as Kl-2). A rough in- 
spection of Fig. 4 suggests that these features might be affected 
by extinction. The CO spiral arm ridges going across the K nu- 
clear bar, create artificially the weat leading arms signature. 
The same applies to the Kl-2 features. Following a similar pro- 
cedure to the one used for the H^ map, we have corrected to 
the first order the K band image by a factor e^*'' . We estimate 
tk from Ak=0.12 A„ (where A„ comes from equation [1]). In 
the corrected map, the Kl-2 features disappear somewhat and 
the leading arms are less visible. 



5. Kinematics 

Fig. 5 shows the CO (1-0) isovelocity contours superposed on 
the Iio map. Although basically a rotating pattern, the veloc- 
ity field is characterized by the presence of strong streaming 
motions associated with the nuclear spiral arms. Deprojected 
onto the galaxy plane, the latter reach 50-60kms~^ within the 
synthesized beam of ~ 2". 

We used the whole data cube to fit both the systemic veloc- 
ity and the dynamical center, as stated in section 2. The global 
velocity field is highly symmetric with respect to Vays and C. 
However, we note that the systemic velocity contour does not 
go exactly across the dynamical center, coincident with the C 
central clump. As it is seen in the p-v major axis diagram 




Fig. 3. I™Q emission contours (levels as in figure 2) arc overlaid on 
Ha (a(top)); b(bottom) I™^ gray scale contours superposed to 
6cm radio continuum emission contours. 



of Fig 6a (see below), this is due to the strongly asymmetric 
wings of the CO emission line towards C. 

In Fig. 6a we show the major axis p-v diagram. The depro- 
jected CO rotation curve (vo6s) derived from Fig.Ga (assuming 
i=32°) implies a steep rise of ISOkms"'^ in less than 2", as- 
sociated with the ultracompact region C, marginally resolved 
by the interferometer beam. Henceforth, the velocity gradient 
there might be higher. The dynamical mass enclosed in the in- 
ner lOOpc is 6.5xlO^M0, the percentage of molecular mass in 
C is ~20%. 

At larger radius (at r=±2 5") we note the existence of 
a hole in the major axis p-v diagram, at a position coincident 
with the interarm crossing, on both sides of the nucleus. There- 
fore we are unable to determine accurately Yobs in this region. 
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Fig. 4. We overlay I^Q(line contours) on the K-band image (gray 
scale) . 



We can either assume a simple solution consisting of interpo- 
lating between the known values of Vohs at r=2" and r=7", or 
just consider Vobs to be undetermined within this region. The 
uncertainties of the inner Yobs arc discussed in section 6.2, in 
particular a comparison with the rotation curve obtained from 
the infrared potential(vcir) will help to clarify what might be 
the "real" trend of rotation curve in the ND. The latter is 
shown to be crucial in the model fitting. Moreover, although 
tentatively, we report the detection of CO emission from two 
symmetric lobes at r=±4" showing velocities highly forbidden 
by the preponderant rotating pattern (see Fig. 6a). Isoveloci- 
ties appear distorted in this region encircling C on both sides 
of the nucleus (sec Fig. 5). This could be the signature of orbit 
crowding in the vicinity of the ilLR, or alternatively out of the 
plane motions. 

Farther away, and coinciding with the crossing of the in- 
ner spiral arms (r= 10-20") the declivity at both sides of the 
p-v diagram might evidence strong streaming motions asso- 
ciated with the arms. However the symmetry of the declivity 
suggest that we sec a real decline of Vrot, reflecting the exis- 
tence of a compact mass distribution in the ND. Disentan- 
gling between the kinematical signatures of the axisymmetric 
and the non-axisymmetric parts of the mass distribution is not 
straightforward. We illustrate below (section 6) how numeri- 
cal simulations help us in separating both contributions. See 
paper I for a discussion of rotation curve at r>20". 

We show in Fig 6b the p-v diagram along the minor axis. 
We assume the eastern side (x>0, eastwards) to be the near 
side. The most remarkable feature is the strength of the ra- 
dial component of the non circular motions at the crossing of 
the spiral arms. The velocity gradient of the streaming, de- 
projected onto the galaxy plane, reaches ~80kms~^ within a 
beam of ~2" (165pc). 

6. Numerical simulations 

Previous numerical simulations have explored the fit of one 
pattern associated with the primary bar and the spiral struc- 



Fig. 5. overlaid on a rotated Iio grey scale image of MlOO, ori- 
ented along the major (y>0, northwards) and minor (x>0, east- 
wards) axes, we plot the first-moment isovelocity contours derived 
from the 1—0 synthesized field. Isovelocities go from v=-100kms~^ 
to lOOkms^^, by steps of 20kms~^, relative to the systemic velocity 
vsi/s=1573kms^^(LSR) (displayed by a thick solid contour). Dashed 
contours stand for v<0. 



turc in the disk of MlOO (G94,S95). The detailed picture of 
the nucleus available now at different wavelengths allows to re- 
fine the model tuning, in particular, to evaluate the feasibility 
of a two independent patterns explanation for the bar within 
a bar morphology and to compare the different modelling ap- 
proaches. 

We perform simulations of the molecular clouds hydrody- 
namics, using a potential as close as possible to reality. This 
numerical code was developpcd by Casoli & Combes (1982) 
and Combes & Gcrin (1985) (sec the detailed description of 
the model in these references) and it was first applied to a real 
galaxy (the interacting system M51) by Garcia-Burillo et al 
(1993). 

We use galaxy plates of MlOO at wavelengths that trace the 
old population, i.e. close to the real mass distribution. The ap- 
proach is different from what K95's work is based on for three 
reasons: first, they use in their simulations a theoretical poten- 
tial only qualitatively similar to MlOO's. When their theoretical 
potential and their rotation curve are compared with the di- 
rectly derived from observations the discrepancies are flagrant. 
Their rotation curve rises much more rapidly (vrot=270kms~^ 
at r=l kpc), and it reaches a much larger plateau value farther 
out {v-rot >350 kms""*^ for r> .5kpc). The position of Lindblad 
resonances arc therefore completely different in their model 
and in the real galaxy. Secondly their simulations use the SPH 
code, which induces much more viscosity in the gas behaviour. 
Henceforth they enhance the action of viscous torques com- 
pared to gravitational torques. Finally, they intend to fit the 
gas response in the whole galaxy disk by a single fast pattern 
(np=70kms~^kpc~^): this is in contradiction with G94,S95 
who fit successfully a slow pattern for the outer disk. In the 
present work two models axe explored: one slow mode or two 
coupled modes (slow-|-fast). 
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Fig. 6. a(top):The ^^CO(l-O) position-velocity diagram along the 
kinematical major axis oriented with a PA=153°. Y oflFset (in arc- 
sec) is measured along the major axis. Gray contours axe linearly 
scaled from 10% to 90% by steps of 10% of the maximum value. 
b(bottom): The equivalent p— v diagram, but now taken along the 
minor axis. X offset (in arcscc) is measured along the minor axis. 
Same equally spaced levels. Dashed lines cross at the dynamical 
center locus on both figures {vsys, x=0, y=0). 

6.1. Gravitational Potential 

We have obtained the potential in the plane of the galaxy from 
the combination of a large-scale red image of the outer disk 
and a small-scale K-band image of the nuclear region. The 
red image is subject to obscuration in the very dusty nucleus, 
while the K-band image is not extended enough to cover the 
whole galaxy disk. The two images were deprojected onto the 
plane of the galaxy, assuming an orientation of PA (position 
angle) = 153° and i (inclination angle)=32°. To combine these 
two images, we performed a Fourier transform at each radius 
of the deprojected images, and obtained up to the 8th order 
Fourier component (in e***). We plotted the radial distribu- 
tion of the coefficients in both images, and found the relative 
scale factor to make these distributions continuous. The con- 
nection between the two images was done at 10" (~800pc). 
Then the combined image was rebuilt from the Fourier coef- 
ficients. Higher order coefficients, from m = 4 to 8 were com- 
pletely negligible inside the nucleus, and their contribution was 
slight (a few % of the m = 2) but visible, in the large-scale disk. 

Through 2D- FFT transform of this recombined image, the 
gravitational potential was obtained in the plane. We run the 
gas simulations in 3D assuming cylindrical symmetry for the 
gravitational forces within the plane (Fa,, Fy). This hypothesis 
is justified considering that the gas thickness (Hgds 50-lOOpc) 



is much smaller than the stellar thickness (Hstars ~lkpc). 
Therefore for the gas itself, we can neglect the variations of 
both Fx and Fy with z, making plausible the local approxi- 
mation of cylindrical symmetry for them. We assume ilstars 
to be constant with radius. The vertical z-forces are derived 
assuming an isothermal stellar disk with a sech'^{z/H) density 
law. The useful grid used for the potential is 512x512, with a 
cell of 1" (« lOOpc), and a maximum radius of « 25 kpc. 

6.2. The rotation curve of the simulations (vdr), compared to 

observations (vobs) 

In a first step, we chose a constant mass-to- light ratio (M/L) to 
compute the rotation curve adopted for the simulations (vcir) 
and compared with the observed one {Yobs)- The disagreement 
was severe, especially for the inner regions where Vcir- >>Vo6s. 
The fit can be obtained adjusting M/L as a function of ra- 
dius. It seems logical to think that M/L is to be lowered in 
the ND because there might be a non- contribution of young 
supergiants within the K band, which have a low mass-to-light 
ratio. Moreover the large amount of dust associated with the 
high neutral gas column densities of the nuclear region, which 
even in the K-band modifies the light profile, also demands a 
modulation of the M/L ratio as a function of radius. 

It is far from straightforward to choose the radial trend of 
the M/L ratio or equivalently the right Vcir- to introduce in the 
simulations code. Imposing Ycir=Yobs, seems unrealistic as it 
implies an anomalous low M/L ratio for the ND (the implica- 
tions of this limit case solution are explored in section 6.3.1). 
Moreover the determination of Yobs is, at places, controversial 
and it is often mixed with strong non-circular motions (see 
section 5.4). 

Still this step is found to be determinant, as the relative 
positions of resonances, and henceforth, the morphology of the 
gas response to the wave, depend strongly, first on the pattern 
speed(s) of the mode(s), and last but not least, on the finally 
adopted Ydr- Indeed the fit of the optimum rotation curve 
for the simulations must be envisaged as an iterative process, 
where both Omegap(s) and Ydr must be varied coherently to 
get a sensible fit. In particular when a two patterns scenario is 
envisaged, we must necessarily find am overlap between the res- 
onances of the fast and slow patterns, otherwise the response 
of the gas departs markedly from the observations. This coin- 
cidence is justified by the theory of modes coupling in spiral 
disks (Tagger et al, 1987; see also Masset and Tagger, 1997). 
These restrictions limit considerably the parameter space of 
the fit when we try the solution of two patterns (see section 
6.3.2). 

Basically we will perform two families of simulation runs: 
a) with one pattern and Ydr ^Yobs at all radii, and hence- 
forth with a strongly lowered M/L ratio in the ND and b)with 
two patterns, dynamically decoupled, for the inner bar and the 
primary outer bar. The latter implies a different tuning of the 
rotation curve in the ND (implying Ydr >Yobs) and conse- 
quently for the M/L ratio (see fig.7). 

6. 3. Best fit solution 

We then simulate the hydrodynamics of the gas in the potential 
computed from the combined K and red images. The gas is 
launched in equilibrium in circular orbits in the axisymmetrical 
part of the potential. 
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Fig. 7. We plot the principal resonance frequencies Q, Q-k/2 and 
Q+k/2, derived from Vdr, i e., the final rotation curve adopted in the 
simulations. We overlay (with a dashed line) the apparent rotation 
curve obtained from the terminal velocities method (voj,s). The loci 
of the resonances are determined by the speed of the two patterns: 
f2y=160kms~^kpc~^ and Qs=23kms~^kpc~^. 
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Fig. 8. Wc display the particle orbits for molecular clouds in the 
region where the bar instability develops, as they are seen, firstly, 
a(top): from the frame rotating at n=23 kms~^kpc~^ , for the solu- 
tion of a single slow pattern and b(bottom): from the frame rotat- 
ing at n=23 kms~^kpc~^ for the best-fit solution of a double pattern 
(Q/=160kms-ikpc-i (from r=0" to 10"), 0^=23 kms-^kpc"! (for 
r>10"). The length of the arrows is proportional to the particle 
speed in the rotating frame. 




Fig. 9. We represent the particle orbits in the ND seen from a frame 
rotating at f2=23kms~^kpc~^, for the slow pattern solution of 6.3.1 
but taking Vdr tuned for the best-fit solution of the double pattern 
of fig 8b 



The decomposition of the stellar potential has been done 
in one axisymmetric part and two non-axisymmetric parts, to 
include the possibility of two nested bars with two different 
pattern speeds. The nuclear bar in the K-band image has a 
radius of about 10'" (~ 1 kpc), and we overlapped the Fourier 
analysis of the images also at 10" . 

The axisymmetric part of the density is first used alone 
to compute the axisymmetric part of the potential; then we 
compute the potential over the whole disk due to the central 
non-axisymmetric part of the density for radii smaller than 
10" , and a third FFT transform is completed for the external 
non-axisymmetric part of the density ( r > 10"). The second 
part of the potential is rotated with f2/, generally at a larger 
angular speed than the main bar, rotating with Qs- Obviously 
when Os=r2/, we recover the one mode solution. 

We have not considered here the self-gravity of the gas, 
since it only slightly modifies the overall morphology, given the 
gas to total mass ratio of NGC 4321 (see S95). We have also 
run a simple collisional code for the gas, without mass spec- 
trum. Clouds interact with each other via inelastic collisions, 
with their radial relative velocity losing 75% of their absolute 
value in the collision. The collisional grid is two-dimensional, 
since the thickness of the gas disk is not larger than the col- 
lisional mean free-path. Clouds are launched with an initial 
radial distribution which is exponential with a scale length of 
7 kpc (70"). 

The non-axisymmetric part (one or two bars), rotating as 
a solid body, is introduced gradually over a time scale of 150 
Myr. Simulations carry on until a steady response of the gas is 
reached; on average the condition is fulfilled after ~ 900 Myr. 

For each numerical run we compute a data cube of the 
simulated galaxy. We use the positions and radial velocities of 
each cloud at the final time of the simulations. The cube is 
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Fig. 10. a(top): Superposition of the simulated gas response (grey 
scale) with the molecular gas distribution seen with the 30in (line 
contours) and the VLA Ihi map (dashed contours) within the opti- 
cal disk of the galaxy, b (bottom): The same but comparing the Ico 
interferometer levels (line contours) to the simulated gas response 
(grey scale). 



synthesized using a convolution beam similar to that of the 
interferometer map, of FWHM~2" for the ND and a 30m-like 
beam of FWHM~13" for the outer disk. 

The gas response is mostly sensitive to the action of grav- 
itational torques induced by the non-axisynimctric potential 
of the disk. Depending on the velocity pattern (s) of the mode 
(s) and the adopted Vcir, molecular clouds accumulate towards 
the loci of the principal resonances of the disk, delineating 
pseudo-rings, ovals and spiral edges. The morphology of these 
substructures and the kinematical signature revealing the pop- 
ulation of different orbit families during the secular evolution 
of the gaseous disk, are to be confronted with the observations 
and henceforth we derive from this detailed comparison the 
optimum fit. 



Fig. 11. Superposition of the particle orbits for the best fit, seen 

from the frame rotating at 17^=160 kms~^kpc~-'^ (arrows) with the 
simulated gas response (grey scale) in the ND of MIOO. Ellipses 
trace the loci of the principal resonances (CORj, oILRy, ilLRy) for 
the jast pattern, in the epicyclic approximation. 

6.3.1. The one pattern solution 

We extend the slow one pattern solution, found by G94 and 
S95 for the outer bar+spiral structure, also to the ND, i.e. the 
two non-axisymmetric parts share the common pattern speed: 

0^=0/^23 kms^^kpc"^ The main difference introduced by 
the model developed in this section is the inclusion of a bar 
component in the potential of the ND, present in the K band 
image. Note that G94 and S95 used a red plate where the 
nuclear bar is not present. Therefore their orbit structure for 
the gas clouds was probably not realistic for the ND. 

The first family of runs assume a fit on Vcir ~Voi,s, fol- 
lowing the approach of S95. In the region where there is no 
information on Vo6a (between r=2-7") we interpolate linearly 
between the observed values. In the epicyclic approximation 
we derive the existence of two ILRs: a) iILR=0.4 kpc and b) 
oILR=2.5kpc (sec Fig.5 of S95). 

The response of the gas for the outer disk (r>30") is similar 
to the best fit solution of S95: the fit of the outer bar-l-spiral 
structure is optimized with f2p=23 kms~^kpc~^. In the ND, 
once a quasi-stationary state is reached after 800 Myr, an el- 
liptical gaseous ring appears between r=3.5kpc and 2kpc, i.e, 
in the vicinity of the oILR. An additional and less conspicuous 
ellipsoidal ring appears between r=0.5 and Ikpc, i.e., in the 
vicinity of the ilLR. There is no trace of gaseous orbits aligned 
along the bar. The particle orbits, as seen by the rotating frame 
at np=23 kms^^kpc^^ are displayed in Fig 8a. 

Tlic formation of rings near the outer and the inner ILRs 
are explained by the long run action of gravitational torques, 
extremely efficient on a cloud-structured medium as it is molec- 
ular gas. The sign of the torques change at the crossing of 
wave resonances (ilLR, oILR, COR and OLR) and also near 
the maximum of the precession frequency Q-it/2. At the lat- 
ter position gas spiral arms, formed by the relative precession 
of the bar driven orbits, change from trailing to leading. The 
angular momentum transfer to gas leads to the formation of 
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rings near the OLR and also associated with the ILRs. A ring 
can be formed at an intermediate position between the inner 
and outer ILR or eventually two rings associated with each 
resonance, as seen here. 

Although the solution is clearly valid for the outer disk, it 
fails completely to explain the observed morphology of the CO 
ND: 

— First the size and the orientation of the rings hardly corre- 
spond with the observed distribution of molecular gas: the 
gas ring formed near the oILR (r-^S Kpc) is much too large 
compared with the size of the CO ND (r~l-1.5 Kpc). 

— We also emphasize the absence of gaseous spiral structure 
in the simulations when they reach steady state. In our 
model, the gas is stopped at the oILR barrier forming a 
ring-like source at a distance of r~3 Kpc, and it cannot go 
inwards. Viscous torques acting on the gas at the oILR ring 
only operate at a much larger time scale and are inefhcient 
to drive the gas inwards, opposed to gravitational torques. 
The morphology is quite similar to the one obtained by 
K95 

— The formation of the inner ring at ilLR is also explained 
by gravitational torques, still the amount of gas clouds 
trapped in the inner 200pc is clearly insufficient and it can- 
not account for the strength of the ultracompact source C. 

If we adopt Ydr directly obtained from a constant M/L 
ratio (vcir »^obs in the ND) the oILR locus is at smaller 
radius (r~l-1.5Kpc), but although the formed ring coincides 
with the size of the ND, no spiral structure is delineated by the 
gas and no ultracompact source is formed at the center (see the 
ND particle orbits in Fig. 9). With the latter Vcir we certainly 
get the gas closer to the center, but not close enough. Moreover 
the response of the gas departs also from the observed one. 

6.3.2. The two patterns solution 

To summarize, the major problem of the one pattern solution is 

its inability to drive the gas to the inner 500 pc and henceforth 
to account for the observed morphology of the CO ND. An 
alternative solution consists of increasing the pattern speed of 
the mode (this is implicitly assumed by K95 who take fip=70 
kms~'^kpc~^). That certainly drives the gas inwards, as the 
principal resonances shrink if we take a larger pattern speed. 
However this solution is incompatible with the outer bar-|-spiral 
structure fit of G94 and S95): a fast mode solution for the 
outer disk can be excluded. In fact, the isodensity contours 
seen in Fig 14 of K95's model are similar to our fig 8a of our 
slow one mode solution. The only difference is the scale due to 
their higher pattern speed. 

The high precession rate of the gas orbits in the ND 
strongly suggests the onset of a fast mode: the maximum of 
Q,-k/2 reaches 200 kms~^kpc~^ and a weighted average preces- 
sion rate (derived from < n-«;/2>= J(f7-re/2)Icodr/ Jlcodr, 
in the inner Ikpc) is close to 150±20kms^"'kpc^^. In the slow 
pattern solution, the rapidly precessing cloud orbits see just 
an ajcisymmetric average of the nuclear bar potential, and fail 
to get in phase with the slow bar forcing. The straightforward 
approach is to try a two pattern solution: one slow pattern 
speed (Qs) for the outer disk and a faster mode (O/) for the 
ND, associated with the K-band bar, and whose speed is to 
be adjusted simultaneously with Vcir- Henceforth we build the 
parameter space of the fit with, first, f2/ and secondly, with 
Vcir (obtained via a tunable M/L ratio). 



We emphasize that the tuning of both cannot be made 
independently: we explore a range of pattern speeds between 
n/=80-180 kms^^kpc^^, and for a fixed fi/ we tune Vdr so as 
to get spatial overlapping of corotation of the fast mode with 
oILR of the slow mode, otherwise the dynamic decoupling of 
modes cannot occur (Tagger 1987; Masset and Tagger, 1997). 
The choice for the range of Qf (80-180) is not arbitrary: in fact 
outside this range there is no possible overlapping of resonances 
and no acceptable solutions are to be expected. We underline 
that the above criterium is a first order approach (strictly valid 
in the epicyclic approximation), and it is taken as a starting 
point: in fact the loci of resonances are not circular- like but 
they can be severely distorted in a barred galaxy(G94). Also 
they have a broad radial extent and therefore by overlapping 
we do not mean mathematical coincidence {^cor=^oilr)- In 
any case, the result of simulations will serve as the conclusive 
test to judge the goodness of the adopted solution. 

Among the explored space of parameters the best agree- 
ment with the observations is found within the interval 
n/=160±20 kms"^kpc"\ Outside this range, i.e. below 
f2/=140kms~^kpc~^, the tuned Vcir fulfilling the overlapping 
condition makes appear anomalous responses in the gas for 
the outer disk, in particular, a strong m=3 pattern outside 
the large-scale bar. Although paper I shows observational ev- 
idences of this pattern (both in the red plate and in the molecu- 
lar gas response), the m=3 mode is weaker than the bar driven 
m=2 pattern. Moreover the rings or structures formed in the 
gas response depart markedly from the observations and no 
steady state is reached in the gaseous disk. 

Fig 8b shows the plot of the particle orbits for the 
best fit evaluated by a standard Chi-squared statistical test: 
fij=160kms~^kpc~^. Figs lOa-b display the corresponding 
synthesized 30m and interferometer-like cubes. The following 
observed features are accounted for by the best fit solution: 

— Formation of a gaseous spiral response in the ND with size 
and pitch angle comparable with the observed in the CO 
mini-spiral. The spiral is formed by precessing xi orbits of 
the fast pattern, mostly inside its corotation (r=14", i.e., 
Rf„, = 1.2kpc). 

— The gas is efficiently driven inwards by the nuclear bar, 
overcoming the oILR barrier of the slow pattern (towards 
r=25", leading to _Rfj£_^=25"): it falls inwards spiraling 
from R'oiLRj crossing R^^ and forming an ultracompact 
source inside the ilLR of the fast pattern (at r=2.5", giv- 
ing Rfjj^ ji=200pc). The strength and unresolved size of the 
central source are in excellent agreement with the observa- 
tions. 

— The goodness of the fit for the outer disk (spiral+large scale 
bar) is maintained. 

— Sakamoto et al, 1995 found that the contribution of the 
nuclear bar to the overall potential is dominant in the inner 
2 kpc region, and that the gaseous spiral structure is mainly 
driven by the imclear bar. The rotation curve fitted in our 
simulations (vcir) shows a declining slope from the center 
up to r=2 kpc. This is an additional argument supporting 
the idea that the nuclear bar is dynamically decoupled from 
the outer bar and that both entities evolve independently. 
The high precession frequency Q-k/2 in the inner 2 kpc 
also suggests the onset of a fast mode for the ND. 

— Fig 11 shows the superposition of particle orbits seen from 
the fast frame with the simulated gas response. Nearly close 
orbits near Rcor shows the characteristic epicyclic shape. 
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Epicycles extend over a broad region (500-800 pc wide in 
radius), and intersect with the ellipsoidal orbits near the 
oILR of the slow pattern. Orbit crowding in the region 
shared by the two families of orbits explain the efficiency 
of the inwards gas transport. Moreover, the structure of or- 
bits sustaining the spiral gaseous wave show that there is no 
net upstream or downstream flow of gas all along the spi- 
ral arms. This, together with the derived high-extinction 
towards the CO maxima, explains why Ha maxima ap- 
pears either upstream or downstream the CO arms show- 
ing no regular distribution, contrarily to the predicted gas 
behaviour in the single slow pattern scenario. 
— There is weak evidence of leading spiral arm structure in 
our simulations. The onset of a leading arm instability ap- 
pears, as theoretically predicted, close to the crossing of the 
ILR of the fast mode (towards r~2-3") but it is washed out 
rapidly to form the ultracompact source C. Note however 
that the leading arms would form at a much smaller radius 
than reported by K95 (they assume a slower pattern for 
the ND). 

7. Summary and Conclusions 

Simulations of the H2 cloud hydrodynamics in the double 
barred system MlOO have shown that the ensemble of observer 
tions (optical, infrared, HI and CO maps) arc best explained 
by a two independent bars scenario. The primary stellar bar 
(of 4.5 kpc radius) and the outer spiral structure share a com- 
mon pattern speed of 513=23 kms^^kpc^^ which places corota- 
tion at R§ofl=8-9 kpc, i.e. beyond the bar end-points though 
well inside the optical disk. Although the nuclear stellar bar is 
mostly aligned with the primary bar (within 20°) it has been 
shown to lead a fast pattern rotating at r2/ = 160kms~^kpc~^, 
having corotation at Rcoh=1-2 kpc radius. Both modes are 
dynamically decoupled and they show overlapping of their ma- 
jor resonances: corotation of the fast mode falls well within the 
ILR region of the slow mode. 

The present model explains the efficient gas transport to- 
wards the nucleus, suggested by the interferometer observa- 
tions, as a consequence of secular evolution driven by the stel- 
lar bar. Molecular gas crosses the ILR region of the slow pat- 
tern, spiraling inwards and forming a trailing spiral structure 
and an ultracompact source encircled by the ILR of the fast 
pattern (Rf^£jj=2.5"). Alternative solutions are unable to ac- 
count for the CO observations. In particular, in the slow pat- 
tern solution gas is stopped at the ILR barrier and forms a nu- 
clear ring outside the ND extent. No central gas condensation 
is formed either. The fast pattern solution proposed by K95 
(Slp=70kms^ kpc^ ) worsens the fit for the outer bar-l-spiral 
structure found by GB94. In addition, two independent meth- 
ods based on the morphology of the residual velocity field for 
the gas (S95) and the identification of spurs in optical pictures 
(e.g. Elmegreen et al 1992) confirm the value of Rcoh reported 
above. 

We conclude that the gas response derived from the CO in- 
terferometer map, and the relation between the different stellar 
(K image) and gaseous tracers of the ND (Ha) are best ex- 
plained by the two pattern model. In particular, it explains 
the high CO concentration in the central part. This gas con- 
centration could be eventually the cause of the nuclear bar 
destruction in this fastly evolving galaxy (see Norman et al 
1996). 
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